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ABSTRACT 

The diffusion of uni-directional magnetic fields by two dimensional turbulent 
flows in a weakly ionized gas is studied. The fields here are orthogonal to the 
plane of fluid motion. This simple model arises in the context of the decay of 
the mean magnetic flux to mass ratio in the interstellar medium. When ions 
are strongly coupled to neutrals, the transport of a large-scale magnetic field 
is driven by both turbulent mixing and nonlinear, ambipolar drift. Using a 
standard homogeneous and Gaussian statistical model for turbulence, we show 
rigorously that a large-scale magnetic field can decay on at most turbulent 
mixing time scales when the field and neutral fiow are strongly coupled. There is 
no enhancement of the decay rate by ambipolar diffusion. These results extend 
the Zeldovich theorem to encompass the regime of two dimensional flows and 
orthogonal magnetic fields, recently considered by Zweibel (2002). The limitation 
of the strong coupling approximation and its implications are discussed. 

Subject headings: ISM: magnetic fields — magnetic fields — turbulence — 
diffusion 
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1. INTRODUCTION 

In the interstellar medium, where the bulk of fluid consists of neutral gas, magnetic 
fields appear to be constantly lost from the fluid — a phenomenon often referred to as 
ambipolar diffusion (Spitzer 1978). It is simply because magnetic fields move with (or are 
tied to) ionized gases while there is slippage between the motion of ionized and neutral gases. 
For parameter values typical of interstellar clouds, the ambipolar drift, however, appears 
to be too slow (roughly by two orders of magnitude) to explain the large dispersion in the 
correlation between density and magnetic field strength (Zweibel 2002). This naturally 
motivates the exploration of other transport processes. In particular, the interaction of 
nonlinear, ambipolar diffusion and turbulent advective mixing is a question of obvious 
relevance. In this Letter, we seek to examine the interplay of these two processes and to 
determine a bound on the turbulent transport of magnetic fields. 

While it is well known that turbulent mixing leads to a rapid diffusion of passive scalar 
fields, this is no longer the case for the diffusion of even a weak magnetic field (far below 
equipartition) in a fully ionized gas (Cattaneo & Vainshtein 1991; Gruzinov & Diamond 
1994), due to the backreaction of the Lorenz force. In a weakly ionized gas, such as in the 
interstellar medium, the problem becomes more complicated since turbulent mixing also 
depends on the collision frequency between ions and neutrals as well as on the strength of 
magnetic fields. It is because in a weakly ionized gas, ions undergo the frictional damping 
due to collisions with neutrals, which effectively reduces the effect of the Lorentz force (or 
Alfven waves). In fact, the previous work by Kim (1997) demonstrated that in two spatial 
dimensions (2D), the diffusion is still reduced below its kinematic value, but the critical 
strength of a large-scale magnetic field above which the diffusion is reduced can be larger 
(by a factor of ^i/j„r, where and r are the ion-neutral collision frequency and correlation 
time of neutrals) as compared to what happens in the case of a fully ionized gas. 
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For simplicity, in this Letter, we consider the mixing of uni-directional magnetic 
fields by two dimensional fiows which are perpendicular to these fields. Note that this 
configuration is different from traditional 2D MHD, where the fields and fiows are coplanar. 
To maintain this geometry, it is necessary to make the strong coupling approximation 
(Spitzer 1978; Shu 1983; Zweibel 1988), by assuming that the drift between ions and 
neutrals is balanced by the Lorentz force on ions due to frequent ion-neutral collisions 
(frictional damping). Our work is the generalization of Zweibel (2002) who considered the 
diffusion of magnetic fields by highly idealized fiows made up of an ensemble of hyperbolic 
stagnation points. Since it is not altogether clear how to relate these fiows to reahstic 
turbulence models, we take a statistical approach here and rigorously derive the diffusion 
rate by assuming a standard scenario of Gaussian and homogeneous turbulence. Note that 
it is possible that this simplified statistical model may fail when the nonlinear, ambipolar 
diffusion is dominant, in which sharp front-like structures are generated (Brandenburg & 
Zweibel 1994). 

Under the strong coupling approximation, magnetic fields are advected passively by 
neutral fiows and diffused by nonlinear ambipolar drift (in addition to the usual Ohmic 
diffusion). Thus, in view of this nonlinear diffusion, one may naively expect that a 
large-scale magnetic field would decay at a rate which is significantly enhanced over the 
turbulent (kinematic) value. We show, however, that it is not the case, because of strong 
fiuctuations (see §3). Specifically, in deriving a generalized Zeldovich theorem (which relates 
the macroscopic quantity (transport) to microscopic dissipation) for a weakly ionized gas 
in 2D motion but containing magnetic fields orthogonal to the plane of motion, we show 
that the fiux transport due to the advection by neutral fiows has an upper bound given by 
kinematic value while the nonlinear diffusion arising from ambipolar drift is insignificant. 
The remainder of the Letter is organized as follows. We present our model in §2 and derive 
a diffusion rate in §3. Section 4 contains the summary and discussion of the Letter. 
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2. MODEL 

In a weakly ionized gas with p ^ Pn'> Ph ^in/^^ni is large, i.e., Vin/vni ^ Pnl Pi'> ^■ 
Here, p, and pi are the density of the bulk of fluid, neutrals, and ions; and and v^i 
are ion-neutral and neutral-ion collision frequencies, respectively. Infrequent neutral-ion 
collisions permit us to prescribe the motion of neutral gas, provided that 

^2 V. 



where v ~ v„ and r are the neutral velocity and correlation time, and va = -B/a/^ttp is 
the Alfven speed with respect to the bulk of the fluid (see Kim 1997). Given a neutral 
velocity, the ion velocity follows simply from the strong coupling approximation as 
'^i — '^n+ [(VxB) X B]/47rpi^'j„. Note here that the strong couphng approximation is vahd 
when 

J^in > i>A, (2) 



where i>A — kB/y/Anpi — kvA\J^in/^ni is the Alfven frequency defined by using the density 
of ion (cf. equation [1]), with k being wave number. When this condition is violated, as 
is likely to be case on small scales due to the high frequency of Alfven waves, the drift 
between ion and neutral motions is no longer balanced by Lorentz force, thereby requiring 
a self consistent treatment of ion dynamics (Kim 1997). 

We consider the mixing of uni-directional magnetic fields (say, B = B{x,y)z) by 
incompressible neutral flows v{x, y) in the x-y plane perpendicular to these fields. Thus, we 
treat the 2D motions of 3D fields. Note that our problem is somewhat similar to that in the 
Goldreich-Sridhar model (Goldreich and Sridhar 1997), which also considers perpendicular 
mixing of anisotropic structures. The evolution equation for the strength of magnetic field 
in this geometry can be written in the following form (see also Zweibel 2002) : 

(at + v-v)g = o;V-[gVQ]. (3) 
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Here, Q = B/ p and a = p/47ri/„j, and the (small) Ohmic diffusivity has been ignored. 
The second term on the left hand side of equation (2) represents the advection by neutral 
flow while the term on the right hand side is the nonlinear diffusion by ambipolar drift. 
The ratio of the effects of these two can be measured by ambipolar Reynolds number 
Rad — vI/Xt — {v'^/v\){Ti'ni)- Here, vl — rjk is the kinematic diffusion rate, and 
At = v\/ Vni is ambipolar drift due to total magnetic fields. Note that Rad can be larger or 
smaller than unity, while still being consistent with equation (1). The question is then what 
the total diffusion rate is in the presence of these two (advection and ambipolar drift) effects. 
Would they act together to significantly enhance the diffusion rate over the kinematic value 
Tjkl To answer this question, we assume Gaussian, homogeneous turbulence and evaluate 
the diffusion rate by using a quasi-linear or 'second order smoothing' closure (for instance, 
see Moffatt 1978) in the next section. The observant reader is no doubt puzzled by the 
fact that equation (3) is two dimensional. The motivation for this simplification is that 
turbulence with k • Bq 7^ will bend magnetic field lines, resulting in its conversion to 
Alfven waves and radiation along Bq. Such fluctuations are intrinsically less effective at 
transporting Q, as much of their energy is expended on bending. However, flute-like eddys, 
with k • Bo = 0, are energetically favored for transport, and also remain correlated with the 
flux tube being transported for a longer time. Thus, the incorporation of eddys with finite 
wavenumber along Bq (i.e., k • Bq 7^ 0) will reduce transport of magnetic fields. 



3. Diffusion rate 

We employ the decomposition of Q into large-scale (Q) and small-scale components 
Q' and assume that flows are on small-scales ((v) = 0). The equations for (Q) and Q' are 
then easily obtained as follows: 

dt{Q) + y -{-vQ') = aV-(g'VQ), (4) 
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dtQ' + v-V{Q) = aV-F, (5) 



where 



F ^ [Q"-{Q'^)]V{Q)+2Q'{Q)V{Q) + {QfVQ' 

+Q'^VQ' - {Q'^VQ') + 2(g) [gw - (gvg')] . 

As can be seen from equation (4), the determination of the diffusion rate requires the 
computation of the flux Fj = {viQ') and the cubic nonhnear term in g. To compute 
the cubic nonhnear term, as well as other nonlinear terms that appear in the following 
analysis, we assume that the statistics of fluctuations are Gaussian and that the turbulence 
is homogeneous. Then, 

(g^Vg) = [{Qf + {Q")]V{Q) . (6) 

On the other hand, the flux Fj is evahiatcd by assuming stationary turbulence. We flrst 
multiply equation (5) by Vi, and then take the average to obtain, 

m) = -T^9,{Q). (7) 



Here, r}k — t{v^)/2 ~ vl is the kinematic diffusion rate; kgff — l/l is the inverse 
of the characteristic scale of fluctuating magnetic flelds; Rad — Vk/^T where 
At = [{'^Ay + WA)]/^ni'i 'T is the correlation time of fluctuating magnetic flelds, 
which is assumed to be comparable to that of neutral velocity. Note that and Rad, now 
deflned in terms of averaged quantities, include both mean and fluctuating components. 
In deriving equation (7), we used {vQ'^) — and (viVj) — 5ij{v^)/2 by assuming an 
isotropic turbulence. Equation (7) states how much flux is transported from large to small 
scales, thereby leading to the decay of {Q). Interestingly, the diffusion rate by advection, 
— Fj/9i(g), has an upper bound given by kinematic diffusion rjk, and becomes smaller as 
Rad decreases. This is because that ambipolar drift 'renormalizes' the correlation time so 
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as to reduce the transport. Thus, it is clear that the kinematic turbulent flux is an upper 
bound on Fj. 

In the case of stationary turbulence, the flux transport is balanced by dissipation on 
small scales as: 

m{Q) = -am' + {Q")mQ'f) = -XTmQ'f) . (8) 

This was obtained by multiplying equation (5) by Q' and then taking the average. Equation 
(8), together with (7), establishes the relation between small and large scale fields as 



-{d,{Q)Y = Mmr)- (9) 



1 + R~AD 

Equation (9) is a generalized Zeldovich theorem for a weakly ionized, strongly coupled 
gas. It gives the relation between mean field and its gradient and relates the macroscopic 
quantity (i.e., flux transport) to microscopic dissipation. Note that the original Zeldovich 
theorem in a fully ionized gas can be recovered by replacing the ambipolar drift by Ohmic 
diffusion r) (Ay — ^ rj) and by taking Rj^j^ = 0, which gives {{diQ'Y) j {di{Q)Y = rjk/r] = Rm 
{Rm is the magnetic Reynolds number). Of course, the situation discussed here is 3D, with 
fields orthogonal to the plane of 2D motion. 



Finally, the (total) diffusion rate of (Q) follows from equations (5) — (7) and (9) as 

Vk 



rjT = 



\ , |V(Q)P 



(10) 



1 + ^AD 

rjT given in this form illustrates the two complementary effects of ambipolar drift on 
diffusion — the first is the reduction of rjT by renormalization of r (the term Rad) and 
the second the enhancement of rjT by nonlinear diffusion (the second term in the square 
brackets). Due to the second effect, it is, in principle, possible that rjT S> rjk- This, however, 
turns out to be very unlikely. To see this, we first estimate rjT in an interesting and more 
relevant case where Rad > 1- 
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When Rad > 1, equation (9) leads to 



~ Rad (> 1) , (11) 



suggesting strong fluctuations. Equation (10) then becomes 

rjT ^Vk- (12) 

Thus, rjT approaches the kinematic value rjk- That is, the diffusion cannot exceed the 
kinematic rate, because of strong fluctuations. Note that in this limit, the field is 

transported as an effectively passive scalars! 

We now look at a less interesting limit Rad < 1 where the effect of turbulence does 
not play an important role. Note that in this limit, the strong coupling approximation 
can easily break down. This can be seen by rewriting the validity condition for the strong 
coupling approximation (equation [2]) as Rad > ^/{T^in)- That is, the smaller Rad (the 
stronger magnetic field), the easier it is to violate the strong coupling approximation. Thus, 
the results obtained in the limit Rad ^ 1 may not be consistent with this approximation. 
With this in mind, we reduce equation (9) to 



((vg')') 



Rad (< 1) , (13) 



|V(Q)P 

for Rad < 1- Equation (13) indicates that to satisfy the stationarity condition (equation 
[9]), {Qy > (Q'^). This follows because the dissipation due to ambipolar drift on small 
scales is too large to be balanced by flux transport. However, as noted previously, the 
strong coupling approximation (which leads to the nonlinear diffusion) is likely to be invalid 
on small scales, especially when magnetic fields are strong. What should happen on small 
scales is propagation of Alfven waves, rather than nonlinear diffusion. Alternatively put, 
when the full dynamics of ions is taken into account, stationary turbulence may still be 
possible even when {Q^ > {Q'"^)- Now, the diffusion rate in this case {Rad < 1) follows 
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from equations (10) and (13) as 

VT ~ Vk/RAD ^ Xt^ {vAf/l^ni ■ (14) 

This is a somewhat expected result in the sense that for small Rad, the diffusion rate 
is set by the ambipolar drift. Our non-trivial result is the observation that in this case, 
though the ambipolar drift due to fluctuating magnetic fields must be negligible in order to 
maintain stationarity. 



4. SUMMARY AND DISCUSSION 

The problem of transport of magnetic fields in the interstellar medium is studied 
by incorporating the effect of turbulence. Specifically, we consider the diffusion of 
uni-directional magnetic fields in the presence of 2D incompressible, turbulent (neutral) 
flows perpendicular to these magnetic fields, embedded in a weakly ionized gas. By 
assuming that the strong coupling approximation is valid on all scales, we compute the 
total diffusion rate of a large-scale magnetic field through a quasi-linear analysis. When 
the turbulence is homogeneous, stationary, and Gaussian, the diffusion rate rjT is found 
to depend on Rad and the level of fluctuations (see equation [10]), with ambipolar drift 
playing two complementary roles (see §3). In particular, when Rad > 1, rjT is shown to be 
at most of order the turbulent rate rjk — vl. In this case, the field is effectively a passive 
scalar. Interestingly, this suggests that even in the strong coupling regime, it is unlikely 
that magnetic fields will diffuse at a rate faster than the simple kinematic value, in spite of 
the nonlinear diffusion operator. In the opposite case {Rad < 1), we demonstrated r]T — \t 
as long as fluctuations are negligible compared to the mean field. Note, however, that this 
limit may not be consistent with the strong coupling approximation. Therefore, our result 
not only confirms the main point of Zweibel (2002) but also puts it on a simple, rigorous 
foundation. 
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The results of this Letter are apphcable to any system with a neutral 
populations and weak magnetic field, such that the strong coupling approximation 
(vj = v„ + [(VxB) X B^/iirpiUin) is valid. However, because of the assumed 
incompressibility of neutral flows, a more uniform loading of magnetic fields discussed 
in the paper is basically due to the diffusion of magnetic fields in a constant density 
background. As the magnetic field diffuses while the density remains constant, system 
progresses toward a state of more uniformly loaded magnetic fields. In the turbulent 
case with Rad > 1, this mass loading uniformization occurs by turbulent cascade of 
magnetic energy (by diffusive mixing) to small scale, where it is eliminated by Ohmic 
dissipation. This uniformization occurs in one large eddy turnover time, as one expects, 
since T„„ = EnTn = EJn/vn = EniU e'/%/^) = {ll^' / e^^Xl/ {1 - a^^)) ~ To. Here, e is 
the energy dissipation rate, and Z„ = aln-i with a ~ 1/2 was used. It is important to realize 
that the breakdown of flux freezing on small scales due to Ohmic diffusion is critical to the 
uniformization of loading. Strictly speaking, our results cannot be directly applied to star 
forming regions, where the compressibility of flows and gravity are crucial. Nevertheless, 
our results imply that the ambipolar drift in turbulent medium can make magnetically 
sub critical clouds supercritical and also that the turbulent mixing can uniformize the 
loading of magnetic field lines on a large eddy turnover time scale. 

In the interstellar medium, the diffusion rate due to ambipolar drift alone is too small 
(by a few orders of magnitude) to explain the observations. Given that Rad based on a 
large-scale magnetic field seems to be larger than unity, turbulent mixing perhaps provides 
a mechanism by which uniformity of the density and the strength of magnetic fields is 
achieved on the eddy turnover time scale. However, ions and neutrals are unlikely to be 
strongly coupled on small scales due to the high frequencies of Alfven waves, thereby 
invalidating the strong coupling approximation. Therefore, the complete answer to the 
problem ultimately requires the self-consistent treatment of ion dynamics. Furthermore, for 
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a better estimate on Rad, some information on the strength of fluctuating magnetic fields 
is needed. For instance, when fluctuations are much stronger than mean fields, Rad based 
on fluctuations may be smaller than the unity, and thus ambipolar drift alone may lead to a 
fast diffusion. Of course, even in this case, the validity of the nonlinear (ambipolar) diffusion 
may become questionable due to the breakdown of the strong coupling approximation. In 
either cases, the relaxation of strong coupling approximation is expected to bring in the 
reduction of the diffusion as magnetic fields are no longer passively advected/distorted 
(Kim 1997). Of course, other effects, such as gravity, the detailed microscale mechanism for 
dissipating magnetic energy, and turbulence intermittency must be considered as well. A 
detailed study of this issue will be addressed in a future work. 

We thank E.G. Zweibel for helpful comments. EK and PHD are supported by the U.S. 
Department of Energy under Grant No. FG03-88ER 53275. 
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